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ABSTRACT 

Context. MHD waves can be responsible for plasma fluctuations and short-term variations of the pulsar emission. 
Aims. We consider the properties of plane and cylindrical MHD waves that can exist in the force-free magnetosphere. 
Methods. Waves are considered by means of a linear analysis of the force-free MHD equations. 

Results. We argue that these particular types of MHD waves can exist in the magnetosphere of pulsars. These waves 
are closely related to the Alfven waves of the standard magnetohydrodynamics modified by the force-free condition and 
non-zero charge density. We derive the dispersion equation fror magnetospheric waves and show that the wave periods 
are likely within the range ~ 10"^ - lO"** s depending on the magnetospheric parameters. 
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1. Introduction 

The magnetospheres of pulsars consist of electron-positron 
plasma with, possibly, some amounts of ions. This plasma 
can affect the radiation produced in the inner region of the 
magnetosphere or at the stellar surface. Therefore, under- 
standing the properties of a magnetosphere is of crucial im- 
portance for the interpretation of observations. The grovif- 
ing observational data on spectra and pulse profiles of iso- 
lated pulsars prompt continued improvement of theoretical 
theoretical models of the pulsar magnetosphere (see, e.g., 
Goodwin et al. 2004, Contopoulos et al. 1999, Komissarov 
2006). Apart from a quasi-static structure, however, various 
non-stationary phenomena in the magnetosphere may af- 
fect the radiation. Electromagnetic waves in pulsar plasma 
have been studied extensively over the past few decades. 
The properties of the low-frequency electromagnetic waves 
are of central importance for understanding the underly- 
ing processes in the formation of the radio spectrum. Early 
studies mainly concentrated on the relativistic plasma flow, 
assuming a cold or only mildy relativistic distribution of 
electrons and positrons in the plasma rest frame (see, e.g., 
[ Melrose 1996 and reference therein). Low frequency waves 
were studied by Arons and Barnard (1986), where many 
of the results of the previous studies were rederived and 
generalized. In all these cases, the plasma was assumed to 
be one-dimensional which is justified for plasma in a strong 
magnetic field. The electrostatic oscillations with a low fre- 
quency have been studied recently by Mofiz et al. (2011), 
who found that the thermal and magnetic pressures can 
generate oscillations that propagate in the azimuthal direc- 
tion near the equator. 

Magnetohydrodynamic waves in the pulsar magneto- 
sphere are studied in less detail. Perhaps the only exception 
are diocotron modes, which are the non-neutral plasma ana- 
log of the Kelvin- Helmholtz modes. These modes have been 



studied extensively in the context of laboratory plasma de- 
vices (see, e.g.. Levy 1965, Davidson 1990, Davidson & 
Felice 1998). However, diocotron modes can occur also in 
charged pulsar magnetospheres. The existence around pul- 
sars of a differentially rotating equatorial disc with non- 
vanishing charge density could trigger a shearing instabil- 
ity of diocotron modes (Petri et al. 2002). In the non- linear 
regime, the diocotron instability can cause diffusion of the 
charged particles across the magnetic field lines outwards 
(Petri et al. 2003). This turbulent charge transport could 
bear on the problem of electric current closure in the pul- 
sar inner magnetosphere. The role of a diocotron instabil- 
ity in causing drifting sub-pulses in radio pulsar emission 
has been considered by Fung et al. (2006). Note that the 
diocotron modes should be substantially suppressed in a 
neighbourhood of the light cylinder where relativistic effects 
become important (Petri 2007). The non-axisymmetric dio- 
cotron instability has been observed in 3D numerical mod- 
elling of the pulsar magnetosphere by Spitkovski & Arons 
(2002). The properties of MHD waves should be very par- 
ticular in the magnetosphere pulsars since the electromag- 
netic energy density in the magnetosphere is greater than 
the kinetic and thermal energy density for typical values 
of the magnetic field. This suggests that over much of the 
magnetosphere, the force-free equation will be a good ap- 
proximation for various MHD phenomena, including waves. 
In this Letter we consider MHD waves that can exist in the 
pulsar magnetosphere. 

2. Basic equations 

Despite uncertainties in the estimate of many parameters, 
plasma in the pulsar magnetosphere is likely collisional and 
the Coulomb mean free path of electrons (and positrons) is 
shorter than the characteristic length scale. Therefore, the 
magnetohydrodynamic description can be justified in this 
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plasma. The partial MHD momentum equation for the elec- 
trons and positrons can be obtained in the standard way - 
multiplying the Boltzmann kinetic equation by the velocity 
and integrating over velocity (sec, e.g., Braginskii 1965). 
The sum of electron and positron momentum equations 
yields the well-known force-free condition if one neglects 
terms proportional electron mass and gas pressure, 

PeE+-jxB = 0, (1) 

c 

where pe = e(np — Ug) is the charge density. Taking the 

difference between electron and positron momentum equa- 
tions, one obtains Ohm"s law in the magnetosphere 



(2) 



where a = e'^UpTe/me is the conductivity and V is the 
plasma velocity. 

Eqs.(l) and (2) can be rewritten as 

PeE* + i J X B = 0, J - ctE* = 0, (3) 
c 

where 

E*=E+-xB, J=j-peV. (4) 

c 

If pe 7^ 0, one can exclude E* from Eqs.(3) and obtain the 
equation containing J alone. 



J H J X B 

Cpe 



■■ 0. 



(5) 



It follows immediately from this equation that J = 0. Then, 
Eqs.(3) yields E* = 0. Hence, we have 



j = PeV , E = X B. 



(6) 



These equations imply that the force- free condition (1) and 
Ohm's law (2) are equivalent to the conditions that the 
magnetic field is frozen-in and the only advective current 
exists in the magnetosphere. Or, in other words the force- 
free condition (1) and Ohm's law (2) are compatible in a 
charged plasma only if the electric current is advective and 
the magnetic field is frozen-in. Note that expression (6) can 
be applied only in the regions where pe ^ 0. 

Eqs.(6) should be complemented by the partial continu- 
ity equations for electrons and positrons 



— — -I- V • (n„V„) 
at 



ot — >- a '-a ' 



(7) 



where r^^'' and Ta are the rates of generation and anni- 
hilation of particles of the sort a and is their partial 

velocity. Since electrons and positrons are generated and 

annihilate in pairs, we have re^' = Fp^' and ■* = Fp \ 
The difference of the positron and electron equations yield 
the charge conservation law, 



(8) 



The sum of partial equations (7) yields the equation for the 
total number density 



dn „ , 
^+V.(nV) 



2r„ 



(9) 



3. Equation for MHD waves in the pulsar 
magnetosphere 

MHD processes in the force-free pulsar magnetosphere are 
governed by Eqs.(6), (8), and (9) complemented by the 
Maxwell equations 



V • E = Aivpe, V X E = - 



15B 

^ X, r. ^ T, l^E 47r. 
V-B = 0, VxB = -— + — j, 

c at c 



j = peV, E = X B, 

c 



(10) 



The force-free MHD phenomena are very particular and 
this point can be illustrated by considering linear MHD 
waves. We assume that the electric and magnetic fields are 
equal to Eq and Bq in the unperturbed magnetosphere. The 
corresponding electric current and charge density are jo and 
Peo, respectively. For the sake of simplicity, we neglect hy- 
drodynamic motions in the magnetosphere and focus on the 
effects caused by electric currents. Linearizing Eqs.(lO), we 
obtain the set of linear equations for waves of a small am- 
plitude. Small perturbations will be indicated by subscript 
1. We here consider the waves with a short wavelength and 
space-time dependence oc exp(zu;i — ik • r), where u) and 
k are the frequency and wave vector, respectively. These 
waves exist if their wavelength A = 27r/fc is short compared 
to the characteristic length scale of the magnetosphere, L. 
For these perturbations Eqs.(lO) takes the form 

k-Ei=47rpei, ckxEi^wBi, k • Bi = 0, 
ck X Bi = 4i7rji - cjEi, ji = Pea^i, 
cEi = -Vi X Bo, w/5,.i=k-ji, 
iivni — ik ■ (rtoVi) = 2Fei. 



(11) 



We study only MHD modes with w 4C cfc because one 
can split electromagnetic and hydrodynamic modes in this 
case. At w ~ cfc, a consideration becomes quite cumber- 
some because electromagnetic and hydromagnetic modes 
are strongly coupled. We focus on understanding the natmc 
and qualitative features of magnetospheric waves, rather 
than the direct relevance for observational implications, and 
a consideration of the simple particular case w -C cA; can 
help in that. For the sake of simplicity, we consider waves 
in the region where Fg « 0. Substituting the frozen-in con- 
dition E = — VxB/c into the equation cV x E = —dTi/dt 
and linearizing the obtained induction equation, we have 



iwBi = zBo(k • Vi) - iVi(k • Bo) - (Vi • V)Bo. 



(12) 



Eliminating ji from two equations in the second line of 
Eqs.(ll), we obtain the following expression for Vi 



Vi 



47rp, 



'eO 



-(ck X Bi -FwEi). 



(13) 



Substituting the expression for Ei obtained from the lin- 
earized frozen-in condition and neglecting terms terms of 
the order of (w/cfc)^, we have 



Vi 



AiTcp, 



-Bo X Vi 



IC 



'eO 



47rp, 



■k X Bi. 



(14) 
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The perturbation of the charge density can be calculated 
from the equation = V • Ei/47r, then 

Pel = 7^[*Bo • (k X Vi) + Vi • (V X Bo)]. (15) 

47rc 

Eliminating Bi from Eqs.(12) and (14) in favour of Vi and 
neglecting again terms ~ (w/cfc)^, we obtain the following 
equation for Vi 

47rpeoVi = -k X [(Vi • V)Bo -iBo(k- Vi) + iVi(k-Bo).(16) 

The scalar production of this equation and vector k yields 
the condition 



(k.Vi) = 0. 



(17) 



This equation implies that the longitudinal waves (with 
k • Vi ^ 0) cannot exist in the force-free magnetosphere. 
Only transverse waves with the velocity perpendicular to 
the wave vector (k • Vi = 0) can propagate in this magne- 
tosphere. We have for transverse waves from Eq.(16) 

47rpeoVi - i-(k • Bo)k x Vi = -k x (Vi • V)Bo. (18) 

Eq.(18) is the basic equation governing linear MHD waves 
in the force-free magnetosphere. Note that transformations 
from Eq.(ll) to Eq.(18) were made taking into account 
terms of the two lowest orders in X/L. The term on the 
r.h.s. of Eq.(18) is, in general, of the order of X/L com- 
pared to the second term on the l.h.s. However, this is not 
the case if k is approximately perpendicular to Bq when 
the term on the r.h.s. becomes dominating. 



4. MHD waves in the pulsar magnetosphere 

Plane waves. Consider initally the case of plane waves with 
(k-'B) ^ X/L where we can meglect the term on the r.h.s. 
of Eq.(18). We have for these waves 



47rpeoVi - i-(k • Bo)k x Vi = 0. 



(19) 



The dispersion relation corresponding to this equation is 



a;2 = 



c^F(k-b) 



(20) 



where flm^ircpeo / Bq and b = Hq/Bq. It is convenient to 
express the characteristic frequency flm in terms of the 
Goldreich- Julian charge density, pcj = ^Bo/2Trc where 
is the angular velocity of a neutron star. Then, we have 
rim = and ^ = Peo/PGJ, and the dispersion equation 
(20) can be rewritten as 



w = ±c(k • b) 



ck 



(21) 



This equation describes the new mode of oscillations that 

can exist in the force-free pulsar magnetosphere. Eq.(21) 
is like the dispersion equation for whistlers in "standard" 
plasma. However, there is a principle difference between the 
considered waves and whistlers since Eq.(21) contains the 
charge density peo = e{npo — rieo), whereas the dispersion 
relation for whistlers is determined by engo- Therefore, the 
magnetospheric waves do not exist if peo = 0, but whistlers 



can propagate in neutral plasma. The frequency of mag- 
netospheric waves is higher because it is generally believed 
that \np—ne\ <C rig in the magnetosphere. Deriving Eq. (21), 
we assumed that oj <C ck. Therefore, the considered modes 
exist if 



> c(k-b). 



(22) 



This condition can be satisfied for plane waves with the 
wave vector almost (but not exactly) perpendicular to the 
magnetic field, for which the scalar production (k • b) is 
small. For example, if the angle between k and Bq is (7r/2 — 
6),d 7r/2, then Eq.(22) is satisfied if 



S < ^n/ck. 



(23) 



Note that generally, magnetospheric waves can exist even 
if (5 ~ 1 but our consideration does not apply to this case. 

Cylindrical waves. Let us assume that the basic mag- 
netic configuration of a neutron star is dipole and consider 
a particular sort of waves that can exist in a neighbourhood 
of the magnetic axis. In this neighbourhood, the field is ap- 
proximately parallel (or antiparallel) to the axis but the 
field component perpendicular to the axis is small. We will 
mimic the magnetic geometry near the axis by a cylindri- 
cal configuration with the magnetic field in the 2;-direction. 
Introducing the cylindrical coordinates (s, ip, z) with the 
unit vectors (e^, e^, 62), we can model the magnetic field as 
Bo = Bo{s)ez- Consider the special case of perturbations 
with the wavevector perpendicular to the magnetic field, 
k = (fcs,fc^,0) where fc^ = m,/s and m is integer. Note 
that even though we used a short wavelength approxima- 
tion deriving Eq.(18), m should not be large for cylindrical 
waves because the cylindrical symmetry of the basic state 
is assumed. For these perturbations, Eq.(18) reads 

47rpeoVi = -k X e^Vu (24) 
w ds 

Taking the radial component of this equation, we obtain 
the dispersion relation for cylindrical waves in the form 



cfcy dBp 
4:npea ds 



(25) 



Cylindrical waves may exist only if m ^ 0. These waves 
propagate around the magnetic axis with the velocity 
{c/'lTrpeo){dBo/ds) and period 



P, 



27r 87r speo 



uj m.c{dBo/ds) ■ '■^^^ 
If we represent the unperturbed charge density as S,Peo, then 



am 



(27) 



where P = 27r/n is the rotation period of a pulsar and 
a = din Bo/dln s. The parameter a depends on the mag- 
netic configuration. We can estimate it assuming that the 
poloidal field is approximately dipole near the axis. The ra- 
dial and polar components of the dipole field in the spher- 
ical coordinates (r, 6, ip) are 

Br=Bp[^y cosO, Be = ^Bp[^y sine, (28) 

where Bp is the polar strength of the magnetic field at 
the neutron star surface and a is the stellar radius (see. 
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e.g., Urpin ct al. 1994). The radial field component is 
much greaterthan the polar one near the magnetic axis 
and, hence, a w d\nBr/d\ns. Taking into account that 

= .s^ + and s <C z in the neighbourhood of the axis, 
wc obtain with accuracy in terms of the lowest order in s/r 
that a = —As^/r^. Then, substituting this estimate into 
Eq.(27), we have 

It turns out that the period of cylindrical waves does not 
depend on the distance from the magnetic axis if s <C r but 

depends on the height above the magnetic pole. Therefore, 
at any given height, pcrtiirbations rotate rigidly around the 
axis with the period _P,„. For example, Pm near the polar 
spot at the surface is approximately 

Pm « 6.3 X 10-6 ^PoToV (30) 

where Po.oi = P/O.Ols. The period is longer at a higher 

height, r > a. 

We have considered the cylindrical waves only for the 
dipole geometry. Note that in real pulsars, the magnetic 
field can depart from a simple dipole geometry and may 
have a very complex topology even with small-scale com- 
ponents (see, e.g., Bonanno et al. 2003, 2005, 2006). The 
mechanism of formation of these complex magnetic struc- 
tures is related to the earliest stage of the neutron star life 
when the star is hydrodynamically unstable. Dynamo ac- 
tion induced by instabilities generates the magnetic field of 
various length scales that range from the stellar radius to 
a very short dissipative scale. These magnetic fields can be 
frozen into the crust that forms when the neutron star cools 
down. Owing to a high conductivity of the crust, the mag- 
netic structures can survive a very long time. For example, 
structures with the length scale > 1 km (magnetic spots on 
the surface) survive during the entire active life-time of pul- 
sars (Urpin & Gil 2004). The complex magnetic topology is 
inevitable in neutron stars and can be very important for 
the magnctospheric structure, particularly, in regions close 
to the star. It is likely, therefore, that cylindrical magnetic 
geometries similar to those considered in this section can 
exist in different regions of the magnetosphere. 

5. Discussion 

Our consideration shows that the particular type of waves 
can exist in the force-free pulsar magnetosphere. These 
waves arc closely related to the Alfvcn waves of the stan- 
dard magnetohydrodynamics but they are substantially 
modified by the force-free condition and non-vanishing 
electric charge density in the pulsar magnetosphere. Like 
Alfven waves, the magnctospheric ones are transverse be- 
cause plasma moves only in the direction perpendicular to 
the wavevector, k • Vi = 0. Our simplified analysis applies 
only to waves with k almost perpendicular to the magnetic 
field, k • B <C kB but, in general, these waves can exist for 
any k. 

The frequency of magnctospheric waves depends on the 
ratio of the true and Goldreich- Julian charge densities, ^. 
The pair plasma in the pulsar magnetosphere is likely cre- 
ated in a two-stage process: primary particles are acceler- 
ated by an electric field parallel to the magnetic field near 



the poles up to extremely high energy and these produce 
a secondary, denser pair plasma via a cascade process (see, 
e.g., Michel 1982). The number density of this secondary 

plasma can exceed the Goldreich- Julian number density, 
ncj, required to maintain a corrotation electric field and, 
hence, the factor ^ can be large. Unfortunately, this factor 
is model-dependent and quite unsertain (see, e.g., Gedalin 
et al. 1998, Usov & Melrose 1995) and can be different in 
different regions of the magnetosphere. Therefore, it is dif- 
ficult to estimate the periods of waves with a high accuracy 
but, likely, they are within a range 10^^ — 10^^ s. 

In our analysis, we neglected the destabilizing effects 
of electric currents and shear (that can exist, for exam- 
ple, because of differential rotation). These effects may lead 
to instability of the magnctospheric waves and exponen- 
tial growth of small perturbations. Perhaps the instabil- 
ity caused by electric currents is more efficient in a very 
strong magnetic field. This sort of instability has received 
considerable attention in the literature on thermonuclear 
fusion. The principal features of the pulsar magnetosphere 
are the force-free condition and non-zero charge density. 
Nevertheless, in general, the current-driven instability can 
generate the magnctospheric waves. Another mechanism for 
generating these waves can be relevant for shear and/or dif- 
ferential rotation. It is well known that Alfven waves can be 
generated in differentially rotating plasma with the mag- 
netic field by the magnctorotational instability (Velikhov 
1959). This mechanism can operate even if the magnetic 
field is strong and the plasma is magnetized (Urpin & 
Riidiger 2005). Likely, this type of mechanism also oper- 
ate in those regions of the pulsar magnetosphere where 
plasma rotates differentially. The generation mechanisms 
of the magnctospheric waves will be considered elsewhere. 
Generation of waves will probably lead to fluctuations of the 
magnetospheric parameters and, hence, the emission of pul- 
sars. Measurements of the time scale of these fluctuations 
can provide information regarding the physical conditions 
in the magnetosphere. 
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